Disentangling Scaling Properties in Anisotropic and Inhomogeneous Turbulence. 
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We address scaling in inhomogeneous and anisotropic tur- 
bulent flows by decomposing structure functions into their ir- 
reducible representation of the SO (3) symmetry group which 
are designated by j, m indices. Employing simulations of 
channel flows with Re^ ~ 70 we demonstrate that different 
components characterized by different j display different scal- 
ing exponents, but for a given j these remain the same at 
different distances from the wall. The j = exponent agrees 
extremely well with high Re measurements of the scaling ex- 
ponents, demonstrating the vitality of the S0(3) decomposi- 
tion. 

Most of the available data analysis and theoretical 
thinking about the universal statistics of the small scale 
structure of turbulence assume the existence of an ide- 
alized model of homogeneous and isotropic flow. In 
fact most realistic flows are neither homogeneous nor 
isotropic. Accordingly, one can analyze the data per- 
taining to such flows in two ways. The traditional one 
has been to disregard the inhomogeneity and anisotropy, 
and proceed with the data analysis assuming that the re- 
sults pertain to the homogeneous and isotropic flow. The 
second, which is advocated in this Letter, is to take the 
anisotropy explicitly into account, to carefully decom- 
pose the relevant statistical objects into their isotropic 
and anisotropic contributions, and assess the degree of 
universality of each component separately. We analyze 
here direct numerical simulations of a channel flow with 
Rca ~ 70 0-|3|. The main conclusion of this Letter is 
that this procedure is unavoidable; in particular it high- 
lights the universality of the scaling exponents of the 
isotropic sector which are presumably those governing 
the universal small scale statistics at very high Reynolds 
numbers. In agreement with recent studies of the this 
subject |^,||we report that different irreducible represen- 
tations of the symmetry group (characterized by indices 
J, m) exhibit scalar functions that scale with apparently 
universal exponents that differ for different j. The ex- 
ponents found at low values of the Reynolds number for 
the j = (isotropic) sector are in excellent agreement 
with high Re results; these exponents are invariant to 
the position in the inhomogeneous flow, leading to rein- 
terpretation of recent findings of position dependence as 
resulting from the intervention of the anisotropic sectors. 
The latter have nonuniversal weights that depend on the 
position in the flow. 

We consider here channel flow simulations on a grid of 
256 points in the stream- wise direction x, and (128 x 128) 



in the other two directions, y, z. We denote by z the di- 
rection perpendicular to the walls and by y the span-wise 
direction in planes parallel to the walls. We employ pe- 
riodic boundary conditions in the span-wise and stream- 
wise directions and no-slip boundary conditions on the 
walls. The Reynolds number based on the Taylor scale 
is Re\ « 70 in the center of the channel (z = 64). The 
simulation is fully symmetric with respect to the central 
plane. The flow correctly develops a mean profile in the 
stream- wise direction which depends only on the distance 
from the wall, Ux{z). The mean profile shows the three 
typical regimes: a laminar linear mean profile inside the 
viscous sublayers, a logarithmic profile for intermediate 
distances and finally a parabolic mean profile in the core 
of the channel. For more details on the averaged quanti- 
ties and on the numerical code the reader is referred to 

Previous analysis of the same data-base [y as well as of 
other DNS |6| and experimental data in anisotropic 
flows found that the scaling properties of energy spec- 
tra, energy co-spectra and of longitudinal structure func- 
tions exhibit strong dependence on the local degree of 
anisotropy. For example, in the authors studied the 
longitudinal structure functions at fixed distances from 
the walls: 



S^p\R,z) =< {vx{x + R,y,z)~Vx{x,y,z)Y 
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where < • • • >z denotes a spatial average on a plane 
at a fixed height z, 1 < z < 64. For this set of ob- 
servables they found that: (i) These structure functions 
did not exhibit clear scaling behavior as a function of 
the distance R. Consequently, one needed to resort to 
Extended-Self-Similarity (ESS) |9| in order to extract a 
set of relative scaling exponents (^{p) = C^(?')/C^(3); (h) 
the relative exponents, C^(p) depended strongly on the 
height z. Moreover, only at the center of the channel and 
very close to the walls the error bars on the relative scal- 
ing exponents extracted by using ESS were small enough 
to claim the very existence of scaling behavior in any 
sense. Similarly, an experimental analysis of a turbulent 
flow behind a cylinder showed a strong dependence 
of the relative scaling exponents on the position behind 
the cylinder for not too big distances from the obsta- 
cle, i.e. where anisotropic effects may still be relevant 
in a wide range of scales. In the following we present 
an interpretation of the variations in the scaling expo- 
nents observed in non-isotropic and non-homogeneous 
flows upon changing the position in which the analysis is 
performed. In particular, we will show that decomposing 
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the statistical objects into their different {j, m) sectors 
rationahzes the findings, i.e. scaHng exponents in given 
(j, m) sector appear quite independent of the spatial lo- 
cation; only the amplitudes of the SO (3) decomposition 
depend strongly on the spatial location. These findings, 
if confirmed by other independent measurements, would 
suggest that the apparent dependence of scaling expo- 
nents for longitudinal structure functions on the location 
in a non-homogeneous flow results from of a superposi- 
tion of power laws each of which is characterized by its 
own universal scaling exponent. The amplitudes of the 
various contributions may depend on the local degree of 
anisotropy and non-homogeneity. 

Our method of analysis is quite simple [^Sj . We start 
by a direct measurement of the longitudinal structure 
functions 



S'-P^ (r^ R) = ( [(M(r^ + R)- u{r'' - R)) ■ R 



(1) 



Note that the two velocity fields are measured at the ex- 
tremes of the diameter of a sphere of radius R centered 
at r'^. Due to the inhomogeneity this function depends 
explicitly on r'^. Due to the anisotropy the function de- 
pends on the orientation of the separation vector 2R as 
well as on its magnitude. The average must be taken 
over different time frames. Typically we have used 160 
time frames for such an average. The time frames are 
separated by about one eddy turn over time. In each 
time frame we also improved the statistics by averag- 
ing over one fourth of the total number of spatial points 
in the plane at fixed z, invoking the homogeneity in the 
span- wise and stream-wise directions, x, y. Thus we have 
finally about 1 x 10^ contributions to each average. 

Having computed S'^p^ {r'^, R) we decompose it into the 
irreducible representations of the SO (3) symmetry group 
according to: 



S^P'>{r'',R) 



Y,S^l{r^,\RmrniR) 

j,m 



(2) 



We expect that when scaling behavior sets in (presum- 
ably at high enough Re) we should find: 



(3) 



In other words, we expect Q the scaling exponent C"' (p) 
to be independent of m. 

The first result that we want to display is that by ap- 
plying the SO (3) decomposition we seem to improve sig- 
nificantly the very existence of scaling behavior. In Fig.l 
we show (i) the log-log plot of the raw structure function 
with p = 4 measured on the central plane with the 
vector R in the streamwise direction, R = Rx, and (ii) 

the fully isotropic sector S^^'^l{r'^,\R\) with the average 
in taken on the sphere centered on the central plane 
= 64. 
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FIG. 1. Log-log plot of the isotropic sector of the 4th order 
structure function Sq^q, vs. R at the center of the channel 
= 64 The data represented by (x) correspond to 

the raw longitudinal structure function, S^'^^r'^ — 64, Rx) 
averaged over the central plane only. The dashed line corre- 
sponds to the intermittent isotropic high-Reynolds numbers 
exponents ^(4) = 1.28. 

It appears that already at this fairly low Reynolds 
number the j = sector shows decent scaling behav- 
ior as a function of R. This is in marked contrast with 
the raw structure function for which no scaling behavior 
is detectable (x symbols in Fig.l). For the raw quan- 
tity the method of "extended-self-similarity" (ESS) is 
unavoidable if one wants to extract any kind of apparent 
scaling exponent. In our analysis we found similar results 
also for higher order structure functions. The scaling be- 
havior is improved dramatically for the components and 
it can be seen even without ESS. Nevertheless we will 
use ESS below for quantitative reading of the exponents 
within every sector. 

The second point we would like to stress is the apparent 
invariance of the scaling exponents belonging to the same 
(j, m) sector with respect to changing the spatial location 
in the fiow. To study this issue quantitatively we resort 
to ESS, and examine the relative scaling of, say, structure 
functions of order n with respect to the structure function 
of order 2 for n = 3,4.... The ESS method is applied in 
each (j, to) sector separately. 

In Fig. 2 we show two typical ESS plots for longitudi- 
nal structure functions of order 4 vs longitudinal struc- 
ture functions of order 2 both at the center z = 64 
and at z = 32 in the sector j = 0. Also, in the inset 
the quality of the scaling can be appreciated by look- 
ing at the logarithmic local slopes of log(S'Q*Q(r'^, \R\)) vs 

log(5'Q o(r'^, |i?|)) as a function of R for the same two 
different central position of the sphere: at the center 
of the channel (r^ = 64) and at one quarter of the to- 
tal channel height (r^ = 32). The two curves give the 
same global relative scaling exponent. The best fits for 
the relative scaling exponents in the sector [j = 0) give 
C^=64(4) ^ ^^=64(4)^^^=64(2) = 1.84 ±0.05 at the center 
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and Co"^^(4) = 1.82 ± 0.04 at r\ = 32. This result is 
remarkable and together with the experimental result of 
ref. 1^ it provides strong evidence for the universality of 
scaling exponent as defined in distinct (j, to) sectors. We 
recall that the accepted value of this relative exponent in 
high-Re experiments is C(4)/C(2) ~ 1.82 ± 0.02 [0. 
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FIG. 2. Isotropic sector: ESS plot of log Sfl vs log Sfl 
at the center, = 64, (bottom curve) and at r^. — 31 (top 
curve), the straight lines are the best fits with slope 1.82. 
Inset: local slope of log Sq'^q , vs. log S^^l as a function of R, for 
two different location in the channel: (+) center of the channel 
= 64, ( X ) one quarter of the total height = 32. The local 
slopes are very close. For comparison we have also plotted 
a horizontal curve corresponding to the accepted anomalous 
high-Reynolds number value, C(4)/(^(2) = 1.82. 

Similarly, but affected from larger error bars, one re- 
covers the same invariance with respect to higher or- 
der moments. For instance we measure ^q^^'*(6) ~ 
Cq^^^(6) = 2.5 ±0.1. As for relative scaling exponents 
of higher j sectors, the scaling is less clean and therefore 
we may only quote qualitative estimates. As an exam- 
ple, for relative scaling exponents of the {j = 2, m = 2) 
and (j = 2, TO = 0) sectors we have CI o^^(4) = 1-1 ± 0.1, 

CifHA) = 1.15 ± 0.1, ci 
1. ±0.1. 



=64(4) = 1.3±0.1, Cir'(4) 
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To underline the quantitative improvement resulting 
from the application of the S0(3) decomposition we show 
in Fig 3 the logarithmic local slopes of the raw structure 
functions S'W(r=,i?x) vs. S^^^r''^,Rx) at = 64 and 
at = 32. Also the logarithmic local slopes of the pro- 
jection on the j = sector at the same two distances 
from the walls are presented. As is evident, the raw 
structure function at the center of the channel and the 
two J = projections are is in good agreement with the 
high- Reynolds numbers estimate C(4)/C(2) = 1.82 while 
a clearly spurious departure is seen for the raw structure 
bmctions at r^ = 32. 



FIG. 3. Logarithmic local slopes of the ESS plot of raw 
structure function of order 4 versus raw structure function of 
order 2 at = 64 (x), at = 32 (★) and of the j = 
projection centered at — 64 (-I-), and at — 32(D). Also 
two horizontal lines corresponding to the high-Reynolds num- 
ber limit, 1.82, and to the K41 non-intermittent value, 2, are 
shown. 

Notice that due to the invariance of Eq.(0) under the 
inversion R —>■ — R all the amplitudes Oj belonging 
to sectors with odd j vanish. Similarly, at the center 
of the channel the symmetry with respect to the center 
Rz —Rz forces all the amplitudes of the components 
with j + TO odd to vanish as well. As a consequence, the 
sector (j = 2,TO = 1) is relevant only when the center of 
mass is not in the central plane. When = 32 we re- 
cover indeed good scaling behavior also for this sector but 
with a relative scaling exponent CIT'^^C-P) slightly larger 
then the relative scaling exponents observed for the other 
j — 2 sectors. This fact, which seems to violate the sup- 
posed foliation in the j index asserted in Eq.(^ is not 
well understood at the moment and it may be correlated 
with the presence of large scale coherent structures (hair- 
pin) oriented at 45° with respect to the walls observed in 
all channel flow simulations | pH . 

Finally, we discuss briefly the determination of the 
scaling exponents associated with higher j sectors. The 
scaling exponent Cj=2(2) was estimated by a number of 
authors on the basis of dimensional analysis [p2|-p^ , and 
the result is Ci=2(2) = 4/3. There is no theoretical knowl- 
edge of the actual value of this exponent with intermit- 
tency corrections. Our direct analysis for j = 2 seems to 
confirm the dimensional expectation, in agreement with 
the previous experimental H finding. 
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L being the outer scale. At such small scales only the 
smallest exponent survives, and this is how the alleged 
universality of the small scales is achieved. 



FIG. 4. Log-log plot of Sfl (★); S'^^,] (x) as functions of R 

at = 32 and of 5*2^2 s-s a function of R (+) at the center, 
= 64. The straight line correspond to the expectation of 
dimensional analysis : (^^=2(2) = 4/3. 
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In figure 4 we show the log-log plot of S"^^ „j(r'^, \R\) vs 
\R\ for (j = 2, TO = 2) at the center of the channel, and 
for [j — 2,m — 2) and {j — 2,ni = 0) at = 32, su- 
perimposed with the straight line with slope 4/3. The 
agreement is quite good. Considering the relatively low 
Reynolds numbers and the fact that the projections on 
the different sectors depend on the non-universal prefac- 
tors aj^m in the decomposition (^) , we think that together 
with the experimental result reported in the present 
finding gives strong support to the view that the scal- 
ing exponents in the j = 2 sector are universal. We are 
not able yet to offer the similar support to the possibility 
that all the scaling exponents in the higher j sectors are 
universal. Such a conclusion calls for additional careful 
analysis of the scaling of higher order structure functions 
and higher j sectors. It is outside the scope of this letter, 
but it is currently under active study. 

In summary, we presented three important results that 
follow from the SO (3) decomposition of the longitudinal 
structure functions measured in channel flow simulations 
10,^: these are (i) The scaling behavior is better defined 
in separated (j, m) sectors. This is in contradistinction 
with the raw longitudinal structure function which fails 
to exhibit any scaling at all. (ii) The isotropic (0, 0) com- 
ponent of the structure functions exhibits a universal 
scaling exponent which is invariant to the spatial loca- 
tion in the flow and the distance from the walls, (iii) 
The j = 2 component exhibits a scaling exponent which 
is compatible with the theoretical expectation and is in 
excellent agreement with the experimental measurement 
1^, indicating universality. 

The picture that emerges is that the higher order sec- 
tors are characterized by scaling exponents that are larger 
than the fundamental exponent in the isotropic sector. If 
this is so, it may explain the decay of anisotropy at small 
scales for high Re flows. In the limit Re^ 00 we ex- 
pect scaling behavior at very small values oi R/L with 
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